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^ This report documents microstructural features of CH3(CH2-0-CH2)mCH3 dissolved 

^ in water and n-hexane for m — 11, 21, and 31. Probability densities for end-to- 

(D end distance, and the associated potential-of-mean-force (pmf), are more revealing 

^ of chain microstructures than are the corresponding results for the radii of gyration. 

^ For water, the pmf identifies three distinct regions: loop-closure, globule, and high- 

^ extension regions. The globule region affirms a water-swollen chain, and is not evident 

Oh 

"i in the n-hexane results. Chain C-atom density profiles from the chain centroid are 

s 

O also different in the water and n-hexanc cases. For n-hexane (but not water), the 

^ density profiles are similar for the different chain lengths when the distances are 

scaled by the observed (Rg^) ■ For water (but not n-hexane) and the smaller chains 
>-> considered, the carbon material exhibits a distinctive enhanced concentration, or 

Ph internal condensation, at the centroid core of the structure. 
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I. INTRODUCTION 



Dispersant materials used in response to oil spill^ undergo a demanding effectiveness 
approval process, and approved dispersants are stockpiled near drilling activity and possi- 
ble spill locations. These factors inhibit continuous improvement of dispersants. Rigorous 
laboratory studies of the molecular-thermodynamic interactions at play in dispersant sys- 
tems are, therefore, particularly welcome. In addition, effective thermophysical modeling for 
these systems should improve the fidelity of laboratory-scale simulations of field-scale and 
ocean-scale oil spill dispersion. A typical prerequisite for development of effective molec- 
ularly specific thermophysical models is computational simulation of these solutions on a 
molecular basis. 

Arguably the most important water-soluble synthetic polymers,^^ (-CH2-CH2-0-)„ chain 
molecules are also soluble in common organic solvents and are intrinsic to the dispersant 
materials in use.^ This report establishes molecular-scale microstructures of PEG oligomers 
in aqueous solutions, information preparatory to molecular-scale theories of the structure 
and function of these species in solution. 

The molecular versatility of PEO/PEG polymers in solution makes them a challenge for 
molecularly specific statistical thermodynamic theory. For example, their conformations 
respond sensitively to the solution environment.'^ Common-sized PEG molecules are dis- 
tinctly helical in n-propanoic, isobutyric, and isopentanoic acid solutions coexisting with 
liquid water.' That helicity emphasizes the contrasting behavior that these chains exhibit 
generic coil structures in aqueous solutions. Though the helices are not evident in bulk liquid 
water, and indeed also not in acetic acid nor in isobutanol nor in n-butanol, helix formation 
does seem to require at least a trace of watei^ 

The solvent plays an intrinsic role in the solution thermodynamics also. The conforma- 
tional sensitivity noted above is associated with a polymer size fractionation between those 
coexisting solutions.'^'Sl Experimental evaluation of a Flory-Huggins interaction parameter for 
PEG /PEG in either water or methanol shows compositional dependence that is substantially 
different in these two cases.^EI 

In many applications,'^ and specifically for dispersants used on oil spills,'^ the PEG/PEG 
chains are decorated with junction, head, or capping groups. The work below anticipates 
those molecular designs, but treats exphcitly only the simplest capping group, (methyl) - 
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CH3. Correlations associated with capping groups serve to focus structural analyses and 
the following work focuses on the radial distributions of capping groups. Alternatively, 
spatial extensions of chain molecules in solution are often characterized by a mean square- 
radius of gyration {Rg^)-, and that information is addressed below for the cases treated. A 
conclusion we draw from the recent exhaustive treatment of longer alkanes in watei'^ is 
that distributions of Rg^ are less informative of molecular structure than spatial correlations 
between chemical groups such as capping groups. The results below should provide further 
clarification of that issue. We note that single-molecule pulling experimentd^'^ are becoming 
a new source of microstructural information, though not yet at the molecular resolution of 
interest in the present work. Those important experiments do not track the radius of gyration 
specifically, but instead manipulate lengths between tethered chemical groups. 

Numerous traditional statistical thermodynamics theories of aqueous PEO/PEG solu- 
tions have been suggested; some representatives were noted in our initial report and we 
do not attempt a review here. An alternative non-traditional procedure for theory develop- 
ment accepts that simulations are necessary for evaluation of molecular theories of liquids, 
that simulations will always be carried-out in any case, and so the simulation work might 
as well contribute to formulation of physical theories. Of course, analyses for those pur- 
poses must express the natural physical concepts under-pinning molecular theories, and 
those analyses must be sufficiently general. We have elsewhere arguecP^ that the recent de- 
veloped molecular quasi-chemical theory does achieve those requirements. That molecular 
quasi-chemical theory has lead over that past decade to resolution of the most basic puzzles 
of hydrophobic effects.'^^'^ It is our intention to establish, for the interesting and impor- 
tant aqueous solutions considered here, some of the groundwork required for fully defensible 
molecular theory. 

As noted already, direct high-resolution molecular simulations are necessary for the theory 
development that is sought. Because these systems are of practical importance, a wide range 
of simulation calculations have been carried-out. Again, we do not attempt a review here, but 
merely note a recent contributiorJ^ that gives access to that substantial available literature. 

The results below were obtained by combination of standard molecular simulation pro- 
cedures appropriate for the sampling requirements of different configurational aspects of 
CH3(CH2-0-CH2)mCH3 chain molecules in water and n-hexane solutions.'^ Further specific 
details are provided in the Appendix. We note here that the combined results support 
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comparisons over a range of molecule extensions, between different chain molecule lengths, 
between different-sized simulation systems, and between water and ra-hexane solvents. 

A curiosity of this simulation work is that with the typical time-scales and length-scales 
it is easy to observe (Fig. [T| molecular-scale break-up of a fluid column to produce droplets, 
and then reformation of a fluid column. We hope in subsequent work to return to these 
observations with further analysis. 




FIG. 1. Formation and then break-up of a fluid column by reconnection across periodic boundary 
conditions. For visual convenience, three images are shown laterally in each frame. In this case the 
droplet/column structure is composed of water with one CIl3(CH2-0-CH2)mCIl3 (m=ll) chain 
molecule (in VDW representation) that here resides near the interface. The majority material in 
the background is n-hexane. 
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II. RESULTS AND DISCUSSION 



From small to large end-to-end lengths, the observed probability density (Fig. [2]) displays 
three distinct behaviors, loop-closure, globule, and high-extension (elastic) regions.!^ The 
loop-closure feature is a realization of a primitive hydrophobic bond,'^ which has been 
extensively studied by simulation but not yet susceptible to experimental measurement due 
to solubility limitations.^- In contrast (Fig. [3|, distributions of Rg^ are less distinctive of 
these microstructural features. 

For liquid water the loop-closure feature is largely independent of oligomer size (Fig. |4]). 
This feature in liquid water is different from what is observed in n-hexane and for non- 
solvated chains (Fig. [s]). These points support the identification of the loop-closure feature 
as an experimentally realizable hydrophobic bond. 

The globule region is clear for the hydrated chains, but not for the case of n-hexane 
solvent (Fig. [6]). Thus in the more benign organic solvent the chains are collapsed relative 
to the more swollen hydrated oligomers. This also supports the interpretion of intrachain 
interactions from the perspective of simple aqueous solution examples. 

The profiles of the density of chain C-atoms, measured relative to the chain centroid, are 
different (Fig. [?]) for the two solvents. For n-hexane, the density profiles are similar for the 
different chain lengths when the distances are scaled by the observed (Rg^)^^'^. This similar- 
ity does not obtain for the case for water. In water and for the smaller chains studied, the 
carbon material exhibits a distinctive enhanced concentration in the middle of the distribu- 
tion, a concentration that is not evident for the case of n-hexane. This interesting internal 
condensation is characterized by a length-scale smaller than [R^)^^"^ . Thus the overall den- 
sity profiles show less similarity on the (Rg^)^^'^ scale, and the internal condensation feature 
is not evident for the largest m=31 chain. 



III. CONCLUSIONS 

The probability density for end-to-end distance (and the associated pmf) for CIl3(CIl2- 
0-CIl2)mCIl3 is more revealing of chain microstructures than are the corresponding results 
for the radii of gyration. For water, the pmf identifies three distinct regions: loop-closure, 
globule, and high-extension regions. In water, the loop-closure feature is similar to a prim- 
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itive hydrophobic bond and is insensitive to chain length. The globule region exposes a 
water-swollen chain, and is not evident in the n-hexane results. Chain C-atom density pro- 
files from the chain centroid are different in the water and n-hexane cases. For n-hexane (but 
not water), the density profiles are similar for the different chain lengths when the distances 

1/2 

are scaled by the observed {Rg^) ■ For water (but not n-hexane) and the smaller chains 
considered, the carbon material exhibits a distinctive, enhanced concentration, or internal 
condensation, at the centroid core of the structure. 

Appendix: Simulation Specifics 

Two different simulation techniques were used for the molecular dynamics calculations. 
The first class of calculations used parallel temperin^^D achieve enhanced sampling of the 
mixing characteristics with water or n-hexane solvent, and of CH3(CH2-0-CH2)mCH3 chain 
conformations. The chain molecules were represented by a generalized amber force field 
(GAFF),^the SPC/E model for water,'231 and optimized potentials for liquid simulations 
(OPLS-AA) were used to describe n-hexane in those simulations.'^ GROMACS 4.5.3 molec- 
ular dynamic simulation packag^^Sl ^^s used for all parallel tempering simulations spanning 
the 256-550K temperature range with 32 replicas (for m = 11, and 21, cases) and 40 repli- 
cas (for m = 31). Long-range electrostatic interactions were treated in standard periodic 
boundary conditions using the particle mesh Ewald method with a cutoff of 0.9nm. The 
Nose-Hoover thermostat maintained the constant temperature and chemical bonds involv- 
ing hydrogen atoms were constrained by the LINGS algorithm. Energy minimization and 
constant pressure density equilibration was performed at 300K to set the constant volume 
conditions, and then production calculations for each replica set were extended to 10 ns. 
Parallel tempering swaps were attempted at a rate of 100/ns, and the temperature grid 
chosen resulted in a success rates of 15-25%. 

The second class of simulations obtained higher spatial resolution in the loop closure 
region using the windowing stratification method.'^^*^ In these simulations a harmonic in- 
teraction between the capping atoms of the GH3(GH2-0-GH2)mGH3 chain was exploited 
to concentrate the sampling of the end-to-end distance to each particular window. For a 
GH3(GH2-0-GH2)mGH3 chain in water or n-hexane, we performed simulations with radial 
displacement coordinate covering the range from 3. 0-10. OA uniformly with 15 windows. Tra- 
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jectories were recorded for lOns/window at the temperature of 297. IK. P{r) was then recon- 
structed over the whole range using the weighted histogram analysis method (WHAM).'2SE3 
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FIG. 2. (upper) Probability density P{r) for end-to-end (methyl-methyl) length for [CH3(CH2-0- 
CH2)mCH3](aq) with m = 21. The normalization of this graph is chosen to suggest the analogy 
with conventional atom-atom radial distribution function in liquids, (lower) Potential of the average 
end-to-end forces showing distinct loop-closure, globule, and high-extension regions. The choice of 
normalization for the upper panel sets the origin of the y-axis of the lower panel, and the same 
convention is followed in both cases. Results here for r < 1.0 nm were are obtained utilizing the 
WHAM procedure. Those high resolution results were matched to overall observation of P(r) from 
molecular simulations with parallel tempering. 
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FIG. 3. For [CH3(CH2-0-CH2)mCH3](aq), distributions of square-radius of gyration displaced 
vertically for visual convenience, extracted utilizing GROMACS analysis tools. Scaling of the 
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abscissa brings the results for different lengths into a similar range : {R^^y = 0.61, 0.84, and 
1.22 nm for ra = 11, 21, and 31, respectively. These values are larger than the ideal [l/^/O) 
proportion of the observed (r'^^'^^'^ 
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FIG. 4. High resolution results for r < 1 nm for several chain lengths from windowing/ WHAM 
calculations. 
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FIG. 5. In the loop-closure region, the [CH3(CH2-0-CH2)mCH3] (m = 21) ohgomer behaves 
differently in n-hexane and water, and the results for n-hexane are similar to those obtained with 
no solvent (ideal gas case). The temperature variation here is small, T = 297 K (water), 302 K 
(ideal), and 297 K (n-hexane). 
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FIG. 6. (upper) For CH3(CH2-0-CH2)mCH3 (aq) chains, at constant density corresponding to 
300K and p = 1 atm: parallel tempering results for different chain lengths and system sizes overlap 
each other on normalized end-to-end distance. Globule and high-extension regions are separated 
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at r ~ (r^) for these chain lengths and system sizes, (lower) Corresponding results for n-hexane 
at constant density (T= 300K and p = 1 atm): Separation of a globule region and a high extension 
region is indistinct here. A Gaussian (parabola) model satisfactorily fits the data in this high 
extension regime. 
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FIG. 7. CH3(CH2-0-CH2)mCH3 chain C-atom density profiles from the chain centroid shows 
higher concentration of polymer C-atoms in water the middle of distribution which is absent for 
n-hexane solvated polymer 
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